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PREFACE 

Part  of  the  RAND  research  program  consists  of  basic 
support  studies  in  mathematics.  In  the  present 
Memorandum,  the  authors  investigate  computer  methods  of 
solving  systems  of  differential  equations  through  the 
so-called  Laplaee-transformation  technique. 

Applications  of  this  technique  are  important  in 
connection  with  the  use  of  invariant  imbedding  in  the 
study  of  time-dependent  transport  processes. 


SUMMARY 


The  usual  analytic  methods  of  inverting  the  Laplace 
transformation — such  as  the  Bromwich— Wagner  Integral  In 
the  complex  plane — are  mostly  Impractical  for  numerical 
work.  This  Memorandum  discusses  a  method  applicable  to 
the  numerical  analysis  of  the  inverse  Laplace  transform 
and  includes  numerical  examples  illustrating  this  method. 
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A  NUMERICAL  INVERSION  QF  THE  LAPLACE  TRANSFORM 
1,  INTRODUCTION 

When  one  is  interested  in  obtaining  numerical 
solutions  of  systems  of  differential  equations,  one 
often  finds  numerical  integration  techniques  to  be  too 
time-consuming,  even  when  employing  high-speed  electronic 
computers.  It  would  then  be  desirable  to  be  able  to 
use  other  techniques  for  solving  such  a  system  of 
differential  equations.  The  Laplace  transformation 
immediately  suggests  itself,  and  its  satisfactory  use  in 
numerical  work  depends  on  the  numerical  analyst's  ability 
to  recover  the  function  f(t),  given  its  Laplace 
transformation  F(s). 

The  usual  analytic  methods  of  inverting  the  Laplace 
transformation — such  as  the  Bromwich-Wagner  integral  of 
F(s)  in  the  complex  plane — are  mostly  impractical  for 
numerical  work.  This  paper  discusses  a  method  applicable 
to  the  numerical  analysis  of  the  Inverse  Laplace  trans¬ 
form  and  Includes  numerical  examples  Illustrating  this 
method.  We  shall  represent  f(t)  by  a  series  of 
orthogonal  functions  and  then  find  the  coefficients  of 
this  expansion  from  the  values  of  F(s)  on  a  set  of 
equidistant  points  on  the  positive  real  line.  For 
examples  of  this  technique  applied  to  other  sets  of 
orthogonal  functions,  see  Papoulis's  original  article  111. 
(Further  examples  can  be  found  in  1 3J - ) 
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2.  PRELIMINARY  DISCUSSION 

Throughout  this  paper  we  shall  assume  that  f (t)  Is 
continuous  on  ( 0,oo )  and  approaches  a  finite  limit  as 
t  -»  oo .  Letting 

(2.1)  F( S')  -  Lap(f)  «  f°°  e~stf(t)dt, 

0 

±. 

we  can  make  the  change  of  variables  r  **  e  ;  then  (2.1) 
can  be  written  as 

(2.2)  F(s)  »  J  1  r0”^  (in  l/r)dr. 

0 

The  function  g(r),  defined  by 
g(r)  »  f (in  1/r), 

is  continuous  on  [  0,1) ,  and  we  can  expand  g(r)  in  a 
series  of  polynomials  that  are  complete  and  orthogonal 
on  (0,1).  Let* 

(2.3)  g(r)  ^  2  c.fl».(r), 

k*0  K  K 

where  {4>k}  Is  the  set  of  these  complete  orthogonal 
polynomials.  If  we  write 


n 


The  notation  of  equation  (2.3)  means  that  2  c.$. 

k»0  K  K 


in  2f2(o,l );  that  Is, 


11m 

n-we 


r 


n 


2  “  g(**) 

lo>0  •  • 


dr  -  0. 
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(2.4)  * Ax)  -  s  a, 

n  j„0  J" 

and  let 

(2.5)  \lmf1  t\(x)i2dx, 

0 


(Of  course#  one  could  pick  an  orthonormal  set  [<X>n),  but 
it  is  more  convenient  to  use  an  unnortnalized  form#  as  we 
shall  see  later. ) ,  then  we  have 


(2.6) 


ckAjc  -  J  \(r)g(r)dr 
0 

k  rl  , 

■  2  a  rJf  ( in  1/r)  dr 


-  aJkP(J  +  1), 


Hierefore,  we  can  determine  the  coefficients  c^ — and 
thus  the  function  f(t)  ■  g(e-t) — by  knowing  P(s)  only 
at  the  Integral  points  1#2# . . . . 

In  general#  by  making  the  change  of  variable 
r  -  e  ,  where  a  is  any  positive  constant,  we  obtain 

(2.7)  F(kff)  -  ~  f  1  r^ffl/o  in  l/r)dr, 

W  vj 

0 


and  by  repeating  the  above  procedure  we  oan  recover  the 
function  f(t)  by  knowing  its  Laplace  transform  only  at 
the  points  a#  2c#....  While  equations  (2.6)  and  (2.7) 
motivate  the  use  of  orthogonal  polynomials  and  damonatrate 


that  f(t)  can  be  recovered  by  knowing  its  Laplace 
transform  at  only  a  sequence  of  equally  spaced  points  on 
the  real  line,  we  shall  not  use  equation  (2.6)  but 
instead  shall  make  use  of  a  more  practical  technique  based 
on  the  orthogonal  set  {<X>n}.  Equation  (2.6)  is  clumsy 
for  numerical  work  since  the  coefficients  a^  get  very 
large  in  magnitude  and  oscillate  in  sign. 


3.  THE  ORTHOGONAL  POLYNOMIALS  - 

The  Legendre  polynomials  Pn(x)  are  complete  and 
orthogonal  on  [—  1,  ll j  that  is,  they  satisfy 

/  1  pJ<*)plc<x>ax  “  !EI  V 

-1 

They  are  easily  computed  by  the  recursive  relation 

(3.1)  (n  +  l)pn+i(x)  “  (2n  +  l)xPn(x)  -  nP^te), 
where  PQ(x)  *=  1,  P-j^x)  ~  x*  ’Hie  ^gendre  polynomials 
are  described  and  tabulated  in  £  -4] . 

We  then  define  the  "shifted  Legendre  polynomials" 
4>n(x)  by 

(3.2)  $n(x)  -  Pn(2x  -  1), 

so  that  (4>n)  is  a  complete  orthogonal  set  on  [0,1], 
and 


(3 


*3)  \mf  Vx)\(y)dx  “  • 


Purthemore,  we  have  the  recursive  relation 


T 


SVriKKW 


— 5“ 

(3.^)  (n  +  l)«n+1(x)  -  (2n  +  1) (2x  -  l)«n(x)  -  nfr^fx), 

where  *0(x)  «  1,  ^(x)  »  2x  ~  1. 

We  now  derive  a  useful  property  of  the  shifted 
Legendre  polynomials,  <J>n.  Define 

(3.5)  an(<T,»)  -  »n(e^t)._8tdt  -  W *„(•■**)). 

0 


for  cr  >  0  and  n  «*  0,1,2,.,..  To  derive  an  expression 
for  Gn  in  closed  form,  we  proceed  as  follows.  Prom 
(2.4),  we  see  that 


(3.6) 


»n(.-<'t)  ■ 


n 


ro-0 


„-mat 


and  hence 


(3.7)  0n(<J,») 


2  a  Lapf.^) 

m-0  ^ 


n 

2 

nwO 


amn 

s-wicr 


y 

s  (s+o)  •  (s+ncr)' 


where  is  a  polynomial  of  degree  £  n. 

Consider  the  case  n  >  1.  (For  n  ■  Ov  we  calculate 
immediately  GQ(a,s)  «  Lap(l)  «  1/s. )  By  t’Ji  change  of 
variables  x  »  e"^,  we  get 

(3.8)  an(o,fca>  -  *  x“^n(x)<bc. 

0 

1c— 1 

For  1  £  k  £  n,  x  is  a  linear  combination  of 
♦o . Vr  further,  because  of  the  orthogonality  of  { 3  on 
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£o,l3 ,  the  integral  in  (3.8)  vanishes.  Thus  we  obtain 

Gn(o,kcO  -  0  for  1  £  k  £  n, 

and  has  zeros  at  s  «  a,2o, . . .  ,ncr.  Hence,  we  have 

Qrt  -  Bn(s  -  c)(s  -  2c)... (s  -  no), 

where  B„  is  some  constant,  and 
n 

(3.9) 

Prom  (3.9),  it  follows  that 

lim  sG  (o,s)  -  (  -l)1^  , 

8-0  n  n 

and  from  (3.7)  we  have 

lim  sG_(a,s)  ■  a^  »  4>  (0) . 
bo  n  011  n 

Wie  recursive  formula  (3.4)  implies  that  $  (0)  »  (— l)n, 
and  therefore  B  «  1,  and  we  have  our  desired  result: 

(3.10)  on(c,s)-i  S  fjg  for  n  -  1,2,..., 

m»l 

00(ff,s)  - 


4.  NUMERICAL  INVERSION  FORMULA 


-7- 

_7t4* 

then  by  letting  r  -  e^v  In  (4.1)  and  taking  the  Laplace 
transformation  of  the  result,  we  would  expect  that 

(4.2)  F(ko)  -  \  c  0  (o,ka), 

-u  n 


since  Gn(ff,kc)  vanishes  for  k  £  n.  A  rigorous  Justi¬ 
fication  of  equation  (4.2)  is  given  in  Sec.  7  of  this 
Memorandum.  Prom  (3.10)  we  have 


(4.3) 


Gn(°'k°)  “  **•  for  n  *  1# •  •  *A-1, 

G'oCtfjkc)  —  . 


Therefore,  from  (4.2),  we  obtain 
(4.4)  F(c)  - 

c k*~l  Cj  i_  i 
,(kB) -«&♦£**&••• 


M for  * 


2,3,.. 


Equation  (4.4)  gives  us  the  following  formula  for 
sequentially  finding  the  coefficients  c  : 


(-4.5) 


'n-1 


OP(o), 

ncP(na) 

n+1 
n=T  * 

n+2  . 
n=? 

n-2 

V  A 

n+(i+i: 

)  n+l 

i  Cj 

1-0  1 

nRiir! 

;  ’  n-i 

naP(na) 

/2n— 1 y 
'  n  ' 

n-2 

-  2 

4 _ r\ 

2n — 1 
*“T“ 


m 


2n~l 


Wierefore,  to  determine  the  first  n  coefficients — 
CQ,C1*  '•"*  *  *cn— 1 — 18  necessary  to  know  the  values  of 


the  Laplace  transform  only  at  the  points 
a, 2a, . . . ,na. 

5«  METHOD  OP  NUMERICAL  INVERSION 

When  inverting  the  Laplace  transform  numerically,  we 
wish  to  represent  f(t)  by  the  sun  X  ) ,  and 

then  determine  the  coefficients  c„  using  formula  (4.5). 
One  should  first  determine  the  optimum  a  to  use. 
Although  in  scans  practical  problems  it  is  often  conven¬ 
ient  to  select  a  before  computing  the  Laplace  transform 
at  a, 2a, . . . ,  we  shall  assume  that  the  Laplace  transform 
we  are  using  is  easy  to  calculate  at  all  points  on  the 
positive  real  line  and  that  we  have  complete  freedom  in 
selecting  a.  We  shall  consider  the  case  in  which  we 
know  that  f(t),  in  addition  to  being  continuous  and 
finite  on  [o,a>] ,  has  the  following  properties: 

f(0)  »  0, 

f(t)  is  monobonic  on  the  positive  real  line. 

In  other  words,  f(t)  behaves  very  much  like 
a(l  -  e^), 

where  a  and  a  are  constants,  a  >  0. 

Sines  we  wish  to  approximate  f(t)  with  as  few 
terms  as  possible  in  the  <&n(e"""ff*)  series,  f(t/ 
should  be  well  represented  as  a  linear  combination  of  a 

—  ,  —  .  ~~ . J' 
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few  terms  of  the  form  e“nat.  Let  us  then  write 

(5.1)  f(t)  -  o(l  - ■«-**).  +  h(t), 

where  h(t)  ^  a2 e”80^  +  a^e""^0^  +  •••. 

We  therefore  wish  to  minimize  the  contribution  of 
h(t).  Writing 

(5.2)  F(s)  “f-ife  +H(s), 

we  shall  minimize  H(s)  for  positive  s.  We  therefore 
desire  that 

«/_\  y  a  a 

pV»)  ^  s  “  s+ff' 

or 

(5.3)  sF(s)(s4cr)  ~  oo. 

Picking  a  (and  a)  such  that  (5-3)  is  a  best  least- 
squares  approximation  for  some  set  si's2' ’ *  * ,sn  of 
positive  real  points ,  we  qbtaln 

(5.4)  a  «  - y, 

NC  -  Ar 

where 

N  N  2 

A-  2  s,F(s,),  B  -  2  s,  F(s..), 

1-1  11  i-1  1  1 

C  —  2  BiVfSi ),  D  -  2  8<^F2(s4). 

1-1  1  1  i-1  1  1 


I 
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(If  we  did  not  have  the  condition  f(0)  «  0,  we  would 
Instead  determine  a,  a,  a»  such  that 

After  using  (5.4)  to  obtain  o,  we  then  can  use  (4.5) 
to  determine  the  coefficients  cn  in  the  representation 

f(t)  ~  Z  cn$n(r),  where  t  «  —  i  In  r. 

We  use  the  recursive  relation  (3.4)  to  evaluate  the  $n 
numerically. 

Since  f  (—  in  r)  is  continuous  on  the  closed 
set  [o,l3,  it  is  known  (see  (53)  that,  al though 
2  °n*n^  nxay  not  converge  for  all  r,  the  second  mean 
of  the  partial  sums  of  the  series  will  always  converge 
pointwise  to  the  function;  that  is,  letting 

(5.5)  s  (t)  -  2  Oj^fe"*15), 

k«0  K  K 


an(1)<fc)  -  m  “„<*)- 


we  have 


(5.6)  11m  o  *2'(t)  -  f(t) 

n-«a> 


for  all  t. 
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For  numerical  work,  c  ^  converges  too  slowly 
for  accurate  determination  of  f .  Rapid  convergence  is 
necessary,  since  it  is  difficult  to  determine  the  cft 
accurately  for  large  n.  A  good  procedure  is  to 
consider  the  weighted  mean, 

2  ksk(t) 

(5.7)  Sn(t)  -  - s 

2  k  KX 

k-1 

which  will  converge  reasonably  rapidly  for  most  practical 
problems  in  which  f(t)  is  well  behaved.  The  next- 
section  of  this  Memorandum  includes  examples  of  numerical 
computations  of  the  inverse  of  the  Laplace  transforma¬ 
tion  using  the  method  described  here. 

6.  NUMERICAL  EXAMPLES 

In  practical  applications,  the  accuracy  of  the 
numerically  obtained  inverse  of  F(s)  is  limited  by 
the  accuracy  with  which  P(s)  is  known.  As  can  be  seen 
from  equation  (4.5),  the  number  of  coefficients  cn 
that  can  be  determined  with  reasonable  accuracy  depends 
on  the  accuracy  with  which  the  quantities  P(no)  are 
known.  The  accuracy  of  the  numerically  obtained  value 
of  f(t)  will  then  depend  on  the  rapidity  of  convergence 
of  the  sequence  6n(t). 

Our  numerical  examples  will  then  be  taken  from 
results  obtained  from  computation  of  the  Inverse  of  the 
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Laplace  transform  of  known  functions  using  the  IBM-7090 
computer.  We  consider  the  following  functions  to 
illustrate  the  results  of  our  method: 


(6.1)  F(s) 


„  ,  .-r/ 2  s 

- ,  for  which  f(t) 

»(«  +  1) 


sin  t  for  t  £  tt/2 
1  for  t  >  tt/2 


i  s  ft  for  t  <  1, 

(6.2)  F(s)  «  -  ■■■■■■'§ — ,  for  which  f (t)  ■«  < 

s*  (l  for  t  >  1. 

We  use  the  least— squares  formula  (5.4)  at  five  pointa-^- 
8^  «  1,  i  **  1,.  ..,5— -to  determine  the  cr  used  in 
equations  (4.5).  The  F(s)  we  used  were  confuted  on 
the  IBM-7090  in  two  ways:  (a)  F(s)  was  calculated  in 
"single  precision"  arithmetic,  which  allows  eight 
significant  decimal  digits  accuracy;  (b)  F(s)  was 
calculated  as  above,  and  a  random  fractional  error  In  the 
Interval  [—10“'’,  +10“'’]  was  Introduced  into  the  result  to 
provide  a  final  F(s)  accurate  to  only  five  significant 
decimal  digits.  In  both  of  these  cases,  we  then 
computed  the  sequence  $n(t)  using  the  full  "single- 
precision"  accuracy  of  the  IBM-7090.  By  computing  the 
inverses  of  several  Laplace  transforms,  we  empirically 
discovered  that  dg  usually  gives  the  best  results  when 
F(s)  is  known  to  within  eight  significant  decimal 
digits,  and  usually  gives  the  best  results  for 
F(s)  known  to  five  significant  digits. 
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The  functions  (6.1)  and  (6.2)  above  test  our  inversion 
method  when  f(t)  has  the  desired  properties  of  mono tonicity 
and  boundedness  but  does  not  resemble  an  exponential., 
Function  (6.2)  has  the  added  difficulty  of  having  a 
l!  sharp  comer."  For  function  (6.1)  we  numerically 
obtain  cr  «  0.87826,  and  for  function  (6.2)  we  get 
0  «  0.84565.  Results  of  our  numerical  inversion  of 
these  functions  are  shown  in  Tables  1  and  2,  and  the 
second  is  exhibited  graphically  in  Fig.  1. 

7.  A  MATHEMATICAL  DISCUSSION 

The  purpose  of  this  section  is  to  justify  equation 
(4.2)  on  a  rigorous  basis.  Again  making  the  change  of 
variables  r  -  e""0^*  let 

g(r)  -  f(-  i  In  r)> 

so  that  (4.1)  becomes 

oo 

g(r)  ~  2  r)* 

n-0  n  n 

Let  <  ,  >  be  the  standard  inner  product  in  J?2[  0,lJ ; 
that  is, 

(7.1)  -  / 1  f1(x)f2(x)dx. 

0 

Hence,  from  the  orthogonality  of  {♦  }>  we  get 

(7.2)  cnAn  -  <4'n,g>. 
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Table  1 


INVERSION  OF  THE  FUNCTION  (6.1) 


t 

f(t) 

~TTa T~ 
5 

~^sr~ 

3 

.  5 . 

0.1 

o. 09983 

0.10018 

0.09950 

0.09982 

0.2 

0.19867 

0.20070 

0.20007 

0.20044 

0.3 

0.29552 

0.29985 

0.29787 

0.29995 

0.4 

0.38942 

0.39668 

0.39198 

0.39701 

0.5 

0.47943 

0.48887 

0.48219 

0.48925 

0.6 

0.56464 

0.57423 

0.56779 

0.57449 

0.7 

0.64422 

0.65122 

0.64746 

0.65129 

0.8 

0.71736 

0. 71908 

0.71966 

0.71896 

0.9 

0.78333 

0.77769 

O.78318 

0.77740 

1.0 

0.84147 

0.82743 

0.83736 

0.82703 

1.1 

0.89121 

0.86896 

0.88218 

0.86853 

1.2 

0.93204 

0.90316  ] 

0.91815 

0.90274 

1.3 

0.96356 

0.93093 

0.94616 

0.93057 

1.4 

0.98545 

0.95319 

0.96732 

O.95293 

1.5 

0.99749 

0.97081 

0.98277 

0.97067 

1.6 

1.00000 

0.98456 

0.99364 

0.98455 

1.7 

1.00000 

0.99514 

1.0010 

0.99526 

1.8 

1.00000 

1.0031 

1.0056 

1.0034 

1.9 

1.00000 

1.0091 

1.0082 

I.OO94 

2.0 

1.00000 

1.0133 

1.0095 

1.0138 

3</o 

1.00000 

1.0175 

1.0041 

1.0178 

4.0 

1.00000 

1.0102 

1.0037 

1.0097 

5.0 

1.00000 

1.0063 

1.0052 

1.0053 

6.0 

1.00000 

1.0045 

1.0062 

1.0033 

7.0 

1.00000 

1.0038 

1.0067 

1.0025 

8.0 

1.00000 

1.0035 

1.0069 

1.0021 

9.0 

1.00000 

1.0034 

1.0070 

1.0020 

10.0 

1.00000 

1.0034 

1.0070 

1.0019 

00 

1.00000 

1.0033 

1.0070 

1.0019 

(a)  Results  obtained  with  F(s)  given 
to  8  significant  digits. 


(b)  Results  obtained  with  F(s)  given 
to  only  5  significant  digits. 


r  Wt  V  ST'  '  V* 
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Table  2 

INVERSION  OP  THE  FUNCTION  F(s) 


0.1 
0.2 
0.3 
0.4 
0.5 
0. 6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 
3.0 
4.0 
5.0 
6.0 
7.0 
8.0 
9.0 

lo.o 

oo 


f(t) 

—JW~ 

5 

0.1 

0. 09503 

0.2 

0.19919 

0.3 

0.31035 

0.4 

0.42236 

0.5 

0.52943 

0.6 

0.62734 

0.7 

0.71354 

0.8 

0.78702 

0.9 

0.84784 

;  1.0 

0.89684 

1.0 

0.93526 

1.0 

0.96456 

1.0 

0.98624 

1.0 

1.0017 

1.0 

1.0122 

1.0 

I.OI89 

1.0 

1.0227 

1.0 

1.0243 

1.0 

1.0244 

1.0 

1.0236 

1.0 

1.0074 

1.0 

1.0058 

1.0 

1.0087 

1.0 

1.0109 

1.0 

1.0120 

1.0 

1.0126 

1.0 

1 . 0128 

1.0 

1.0129 

1.0 

1.0130 

8  significant  digits. 


g8 _ 

0.09971  0.0 
0.20221  0.1 
0.30070  0.3 
0.40346  0.4 
0.51122  0.5 
0.61827  0.6 
0.71736  0.7 
0.80294  0.7 
0.87224  0.8 
0.92498  0.8 
0.96268  0.9 
0.98779  0.9 
1.0031  0.9 
1.0112  1.0 
1.0144  1.0 
1.0146  1.0 
1.0130  1.0 
1.0107  1.0 
1.0082  1.0 
1.0060  1.0 
1.0041  1,0 

1.0060  1.0 
1.0034  1.0 
1.0009  1.0 
0.99959  1.0 
0.998 95  1.0 
0.99866  1.0 
0.99854  1.0; 
0.99844  1.0 


s;  given  to 


0.09527 

0.19932 

0.31023 

0.42210 

0.52918 

0.62720 

0.71357 

0.78721 

0.84815 

0.89722 

0.93565 

0.96492 

0.98652 

1.0019 

1.0123 

1.0188 

1.0224 

1.0240 

1.0240 

1.0231 

1.0070 

1.0062 

1.009& 

1.0123 

1.0136 

1.0141 

1.0144 

1.0145 

1.0146 


,  (6)  Results  obtained  with  F(s)  given  to 

only  5  significant  digits. 


0  0.2  0.4  0.6  0.8  1.0  1.2  1.4  1.6 

t 


Fig.  1 


Inversion  of  the  function  F(s) 
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recalling  from  equation  (3.3)  that 

•J 

A  »  ^4>  •  $  ?*  x"‘  -  1  » 

n  s  n*  2n  +  1  * 

Equation  (3.8)  takes  the  algebraic  form 

(7.3)  tfGn(c,lca)  -  <xk”1,4>n>, 
and  (2.7)  takes  the  form 

(7.4)  crp(ko)  -  <xk"1,g>. 

Since  {4>n}  is  complete  and  orthogonal  in  ^2lo,ll, 
and  since 


Gn(cf*ka)  »  0  for  n  >  k, 

equation  (7.3)  implies  that 

.  ,  k-1  cr(L  .(<r,kcr) 

(7.5)  xfc_1.  Z  -4— -  tlx) 
n«0  An  n 


for  k  •«  1,2,.... 

Therefore,  we  have 


.k-1  aG  (o,ka) 

P(lca)  -  ±  <  z  -Ar - ®„,g> 

n**0  %.  n 


k-1  Gn(a,ka) 


k-1 


n»0 


4—  <4n'g>  ■  nfQ  "AtW). 


which  is  our  desired  result. 

Equations  (7.3)  and  (7.5)  yield  an  interesting 
result.  We  let 


aJS**.' 
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Gn(k)  -  an(l,k)  »  aGn(a,ka), 


and  (7*5)  becomes 


k 

2 

n«0 


°n<k 


+  i) 


n 


Vx)  • 


We  therefore  have  obtained  the  formula  for  the  expansion 

If 

of  x  in  terms  of  the  shifted  Legendre  polynomials: 


(7.6)  x^1  -  |  <t0(x)  +  |  ^(x) 


.  5  k  —  1  k  —  2  A  , 

+  E  *  k"'+"T  *  ITT"?  ®2  'x'  +  * 


2n  +  1  k  —  1 

— ic  *  ftt 


k  - 

FT 


2 

F 


♦  •  • 


k  -  n 
FTn 


<X> 


n 


(x) 


+ 


.  2k  —  1  k  —  1  k  —  2 
+ - F~  *  FTT  *  FT? 


SF^hr  \-l(x) 


The  functions  Gn  are  related  in  the  following  way  to 
the  coefficients  Introduced  in  (2.4).  Letting 

G  (m  +  1) 

h  MS  I  .*f..  ■■■■■■.  - 

nm  * 


so  that  the  quantities  are  the  coefficients  in 

(7.6),  we  have  the  following  two  equations  of  the 
"basis  change"  between  (xn)  and  On): 

_  m  m  „ 

t7-7)  4n.*J0%T 


We  therefore  can  use  the  functions  Gn(m)  to  determine 
the  coefficients  a^,  and  conversely. 
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